Doping dependent isotope effects of the quasi-lD electron-phonon system: 
comparison with the high-temperature superconductors 
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The weak-coupling quantum phase diagrams of the one-dimensional (ID) Holstein-Hubbard and 
Peierls-Hubbard models are computed near half-filling, using a multi-step renormalization group 
technique. If strong enough, the electron-phonon interaction induces a spin gap. The spin gap, 
which determines the superconducting pairing energy, depends strongly on the band filling and 
decreases monotonically as the system is doped away from half-filling. However, the superconducting 
susceptibility exhibits a different doping dependence; it can vary non-monotonically with doping 
and exhibit a maximum at an "optimal" value of the doping. For a quasi-lD array of weakly 
coupled, fluctuating ID chains, the superconducting transition temperature T c exhibits a similar 
non-monotonic doping dependence. The effect of changing the ion mass (isotope effect) on T c is 
found to be largest near half-filling and to decrease rapidly upon doping away from half-filling. The 
isotope effect on the spin gap is the opposite sign as the isotope effect on T c . We discuss qualitative 
similarities between these results and properties of the high-temperature superconductors. 



Recent experiments in the high-temperature supercon- 
ductors suggest the presence of a strong, ubiquitous, yet 
unconventional electron-phonon (el-ph) interaction |jj.y| 
It has been known for some time that these materials ex- 
hibit large and strongly doping dependent oxygen isotope 
effects 0, S E] that cannot be described with the con- 
ventional BCS theory used for ordinary metals. In this 
Letter, we show that the quasi- ID electron gas coupled 
to phonons exhibits highly unconventional doping depen- 
dent isotope effects that are qualitatively similar to those 
observed in the high-temperature superconductors. 

The best understood non-Fermi liquid is the spin- 
charge separated one-dimensional electron gas (1DEG) 
0. The properties of the 1DEG coupled to phonons are 
dramatically different from a conventional metal coupled 
to phonons |3, la, 111 [10( ■ Unlike in a Fermi liquid, in ID 
the el-ph interaction is strongly renormalized, and the 
renormalization is strongly affected by direct electron- 
electron (el-el) interactions. Due to these renormaliza- 
tion effects, a weak, retarded el-ph interaction is capable 
of inducing a spin gap and causing a divergent supercon- 
ducting susceptibility, even when the el-el repulsion is the 
dominant microscopic interaction Unlike the case in 
3D, in ID this can occur without a large amount of re- 
tardation. In contrast to BCS theory, the pairing energy 
and superconducting susceptibility are very sensitive to 
the band filling. 

In the cuprates, both the superconducting transition 
temperature T c and the isotope effect exponent ax c = 
—d\nT c /d\nM, which describes changes in T c induced 
by changes in the oxygen mass M, exhibit highly uncon- 
ventional (i.e. non-BCS) doping dependencies. In BCS 
theory, T c is only weakly dependent on the carrier con- 
centration, and aT c has the universal value of 1/2. In the 
cuprates, T c exhibits a maximum as a function of dop- 
ing, and the isotope effect is not universal-it is strongly 
doping and somewhat material dependent. For the un- 



derdoped cuprates, «t c ~ 1 (indicating that, at least in 
this region, phonons play an important role in the super- 
conductivity). As the doping increases, ar c decreases, 
usually dropping below O.f near optimal doping Q. The 
origin of this behavior remains a mystery of high-T c su- 
perconductivity that any successful microscopic theory of 
pairing should explain. Below we provide a microscopic 
theory that is capable of rationalizing it. 

In this Letter, we compute ar c for a quasi-lDEG cou- 
pled to phonons, under the assumption that charge den- 
sity wave (CDW) order is dephasedby spatial or dynamic 
fluctuations of the ID chains [111 Il2fl . For many choices 
of the parameters, «t c is larger than the BCS value at 
small dopings, then drops below the BCS value as the 
doping is increased - the same behavior observed in the 
cuprate superconductors. We show that the quasi-lDEG 
coupled to phonons displays a strongly doping dependent 
T c that can exhibit a maximum as a function of doping. 
This behavior occurs despite the fact that the pairing en- 
ergy, determined by the spin gap A s , is a monotonically 
decreasing function of increasing doping. We also com- 
pute the isotope exponent qa s = — d\n A s /d\nM. We 
find aA, < 0, the same sign as the isotope effect observed 
on the pseudogap temperature in the cuprates 0. 

The technique we employ is the multi-step renormal- 
ization group (MSRG) method described in detail in a 
previous paper [{| . This method treats el-el and el-ph in- 
teractions on equal footing and properly treats the quan- 
tum phonon fluctuations. In it, we start with a micro- 
scopic el-ph model and integrate out high energy degrees 
of freedom, via an RG transformation. This is done in 
multiple steps, as elaborated below. At low energies one 
obtains an effective field theory that is the same as the 
ordinary 1DEG, except for phonon induced renormal- 
izations of the el-el interactions and bandwidth. The 
accuracy of this analytic technique in computing weak- 
coupling phase diagrams was demonstrated in Ref. by 
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comparison with exact numerical results. We shall ap- 
ply it to two models of interacting, spinful ID electrons 
coupled to phonons: 

The ID Peierls-Hubbard (Pei-Hub) Hamiltonian is 

^Poi-Hub = - ^2 [t - l(u i+ i - lij)] (ct CT c i+ i )(r + H.C.) 



E 



Pi ^ t \ 2 

2M 2 y + ' 



+ U 2jni,i-ni4 . 



In this model, acoustic phonons couple to electrons by 
modifying the bare hopping matrix element t by the el- 
ph coupling strength 7 times the relative displacements 
u i+ i — m of two neighboring ions 0. The last term is 
the Hubbard interaction. For this model, we shall ap- 
proximate the phonon dispersion by its value at the zone 
boundary of 2\J k/M = ljq, since the el-ph interaction 
vanishes at zero momentum transfer. 

The ID Holstein-Hubbard (Hoi-Hub) Hamiltonian is 



Pi K 2 
-^- L - -I a 

2M 2 qi 



^Hoi-Hub = -iy^(c| .cj + i ](T + h.c.) + 

i,cr i 

+ gy / 2Mujo S^ifii + U22 n iA n i,l ■ 

i i 

Here dispersionless optical phonons with coordinate qi 
and frequency ljo = \J K/M couple to the electron den- 
sity rii = cJ ff Ci )(T with el-ph coupling strength g |l4|. 
It is convenient to define the dimensionless quantities 

A Pci = 2iV (7sin/c F ) 2 /K, A Ho i = N g 2 /uj , 
U = U/(ttv f ) , S = Mji/wq)/ HE f /u; ) , 

where Nq = 2/(ttvf). As usual, Vf, k Fl and Ep are 
the Fermi velocity, momentum, and energy, respectively. 
(We have set Ti and the lattice parameter equal to unity.) 
H is the chemical potential measured with respect to its 
value at half-filling. In this Letter, we study the range 
uj < fi < Ep (0 < S < 1). The doping concentration 
relative to half-filling is given by x ~ Nou>o(Ep/uJo) S ■ 
Since the MSRG method is perturbative, it is only quan- 
titatively accurate for X V ,U -C 1 , but is believed to be 
qualitatively accurate for A„,t/ < 1 9J. (The subscript 
v stands for Pei or Hoi.) 

Fig. ^presents the A„ — 5 phase diagrams of the above 
models, computed with MSRG for several fixed values of 
U. The phase boundaries separate regions where vari- 
ous types of order have divergent susceptibilities in the 
low temperature limit. The susceptibility that diverges 
most strongly (i.e. dominates) is shown without paren- 
thesis; parenthesis indicate a susceptibility that diverges 
less strongly. The charge sector is gapless everywhere in 
the phase diagrams. Above the thick solid line, the sys- 
tem is spin -gap ped and described as a Luther-Emery liq- 
uid (LEL) |l5|. Below this line, it is a gapless, quantum- 
critical Luttinger liquid (LL) [g. 
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FIG. 1: T — phase diagrams for the Hoi-Hub model (pan- 
els a and b), and Pei-Hub model (panels c and d). (a) and 
(c) are for U = 0.1; (b) and (d) are for U = 0.4. For all 
diagrams, Ef/ujo = 5. Parenthesis indicate a sub-dominant 
susceptibility. SDW stands for 2&f spin density wave, CDW 
stands for 2k F charge density wave, SS stands for singlet su- 
perconductivity, and 4&f stands for 4k F charge density wave. 



The thick solid line in Fig. ^is determined by integrat- 
ing out degrees of freedom between Ep and ojq, and then 
requiring that the total effective backward-scattering in- 
teraction g\ ot (ujo) — U /(1 + UIq) — Ai(cj ) is zero. Here 
In = \n(Ep/ujo) and Ai(wo) > is the effective strength of 
the backward scattering (momentum transfer near 2kp) 
portion of the el-ph interaction. Below the thick line in 
Fig- E g\ ot (u!o) > and the RG flows carry the effective 
g\ ot to zero at low energies, signifying the stability of the 
LL fixed point. Above the thick line, g\ ot (ojo) < and 
the RG flows carry g\ ot to minus infinity at low energies, 
indicating the existence of a spin gap. 

Ai(wo) is determined in two steps: First, one integrates 
from Ep to fi using the RG flow equations that govern 
half-filled systems, resulting in an effective Ai of [9| 



or 



A 



Hoi 



1 - cUl 



l-X no iX/Uj V (1 + cUlo) 3 



At 



1 



(1) 



(2) 



l-\ Pci Y/Uj y [1 - (c[/Z ) 2 ] 3 
for the Hoi-Hub and Pei-Hub models respectively, where 



X 



4 1- yJ(l-cUl )/(l + cUl ) - 2arcsin(cf77 ), 
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Y = 2cUl a /^\ - {cUl ) 2 , and c = 1 - 6. Next, Ai(/x) is 
used as the initial value to integrate from fj, to u>o, employ- 
ing the RG flow equations that govern incommensurate 
systems, resulting in || 



Peierls-Hubbard 



Ai(cJo) = 




(3) 



for either model, where Z = J Q <fuyexp(u)/(l + u) 3 . 
The condition g\ ot (ujo) — then determines the following 
critical values for the microscopic el-ph couplings: 



A 



Gap 
Hoi 



\ Gap 

A Poi - 



where we defined Si = (1 

^2^3/2 



U{[(1 + Ul )S 3 + Z] /Si + X}~ , (4) 

u{[(i + m )s 3 + z] i s 2 + y}~\ (5) 

cUl ) 3 / 2 (l-cUl )- 1 / 2 ,S2 = 

1 - {cUl fY'\ and S 3 = e sm °/ 2 (l + SUl Q )- 3/2 . The 
system is a spin-gapped LEL for X v > A^ ap . 

In the LEL phase, the portion of the singlet 
superconductivity (SS) and 2kp CDW susceptibility 
that is potentially strongly divergent as T — > is 
given by X ss = (m) F )- 1 (A,/E F )(T/E F ) 1 / K ?- 2 and 
Xcdw = (irv F y 1 (A s /E F )(T/E F ) K f- 2 respectively, 
where the spin gap is A s = cj eexp[— l/|^ ot (o;o)|] 
@, and the effective charge Luttinger exponent af- 
ter integrating out states between Ep and u>q is 
Kf 



2g\ ot 



■#'W]/[2 

Here^ ot M = < ? ^)-A 1 (w )+2A 2 and g\ ot = U/2-X 2 , 
where the forward scattering el-ph interaction A2 is given 
by A2 = And for the Hoi-Hub model and X2 — for the 
Pei-Hub model. The contribution of the Hubbard inter- 
action to gl ot (uj ), given by gf(fj.) = -U/(l - cUl ), is 
obtained by integrating out states between Ep and /z, 
since this contribution is unrenormalized below /1. Note 
that Kf is the effective value at low energies because 
g* ot is not further renormalized below loq 16]. 

The SS susceptibility is the dominant one if Kf > 1. 
The thin solid line in Fig. ^1 is the critical line deter- 
mined by — 1, given by 



A Hol - U 



B±y/B 2 - SiAC/2 



(6) 



where A = (SxX + Z)~ l , B = [{2Si - S 3 )A + C7]/4, and 
C = (1 — cIIIq)^ 1 . xss is dominant if the two conditions 
Si AC < 2B 2 and A^,f < A Hoi < A | S ol + are met. The 
thin solid line is absent in Fig. ^ because Si AC > 2B 2 
everywhere. This line is never present in the Pei-Hub 
model with U > 0, because then < 1 always. 

For 1/2 < Kf < 1, xss is still divergent as T -> 0, 
but, for a single chain 1DEG, xcdw is more strongly 
divergent. The dashed lines in Fig. ^are determined by 
K° — 1/2, which leads to the following critical values 
for the microscopic el-ph interactions 
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FIG. 2: Doping dependence of the dimensionless singlet 
superconducting susceptibility xss (thick solid lines), CDW 
susceptibility xcdw (dashed lines), and spin gap (thin solid 
lines), for the Pei-Hub model with U = 0.4, Ef/loq — 5, and 
various values of Ap e i (labeled in plot), xss and xcdw were 
computed for T/ujq =0.1. 



where D = A[S X {A - 5EX) - 5(EZ + S 3 )}/8 and E = 
(6/U + 3)/5 - C. If X v > A° DW , then Kf < 1/2, which 
means that SS is not divergent. 

Examining Fig. QJi, we see that for moderate values of 
Ap i, for example near Ap c i ~ 0.2, xss is not divergent 
near 5 = 0, where Kf < 1/2, nor is it divergent near 
5=1, where A s = 0. However, xss is divergent for a 
certain range of moderate S. Therefore, in these cases, at 
fixed T <C A s , xss exhibits a peak as a function of 6 at 
an intermediate value of S. This peak is shown explicitly 
in Figs. [2] and |3 where we plot xss = t^fXss (thick 
solid line) versus 6 at T/ujq = 0.1, for representative pa- 
rameters. The CDW susceptibility xcdw = ttwf xcdw 
(dashed lines) does not exhibit such a peak. 

In Figs. [5]and|31we also plot A s /u>o (thin solid lines), 
which shows that at low dopings, xss increases with in- 
creasing doping, despite the fact that the superconduct- 
ing pairing strength A s decreases! The reason for this 
discrepancy is the different doping dependencies of the 
effective interactions in the charge and spin channels, 
which determine Kf and A s respectively. It is worth 
mentioning that in the cuprates, the superconducting gap 
also decreases with increasing doping, which in the under- 
doped region occurs at the same time that T c increases! 

We now consider an array of weakly coupled quasi-lD 
chains with dephased CDW, and treat the interchain cou- 
pling J on a mean-field level, which means that T r is de- 
termined by the temperature at which 2 Jxss = 1 |l7lll8j . 
(The numerical prefactor 2 is determined by the number 
of nearest neighbor chains.) In this case, T c exhibits a 
peak at the same 5 where xss is peaked (assuming J is 
doping independent). The isotope effect exponent ap c is 
readily computed |l9| . and is shown versus 6 in Fig. 
for various values of J = J/(itvf). At low dopings, «t c 
is larger than the BCS value, then drops below 1/2 as 
5 is increased. Fig. 0] also shows q;a s , which is weakly 
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FIG. 3: Doping dependence of xss (thick solid lines), xcdw 
(dashed lines), and A s /u>o (thin solid lines), for the Hoi-Hub 
model (panels a and b), and the Pei-Hub model (panels c and 
d), both with U — 0.1 and E f /lu = 5. Xss and Xcdw were 
computed for T/luq = 0.1. The values of Ah i and Ap c i are 
labeled above each plot. 
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FIG. 4: Dependence of the isotope effect exponents qt c 
and «a s on the doping parameter 5 and interchain coupling 
strength J for the Hoi-Hub model with Ah i = 0.275 (panel 
a) and Pei-Hub model with Ap c j = 0.25 (panel b). For both 
panels, U = 0.1 and Ep/uo = 5. qa s is independent of J. 



doping dependent and negative. Note that there exists 
a range of S for which \<xt c /cxa 3 \ *C 1. Similarly, near 
optimal doping in the cuprates, \aT c /ctT>\ <§C 1, where 
aT , — —d\nT*/d\nM < 0, and T* is the pseudogap 
temperature 0]. 

In some models of high-temperature superconduc- 
tivity based on stripes [2(| , the concentration of holes on 
a stripe remains fixed when the doping in the Cu-0 plane 
changes, but the spacing between the stripes changes. In 
that case, as the doping increases, the parameter 5 re- 
mains fixed, but J increases due to the decreased spac- 



ing between stripes. Then Fig. 0] predicts that ax c again 
decreases with increasing doping. In such a model, in 
the underdoped region, where the stripes are far apart 
and represent well defined quasi- ID chains, increasing 
the doping increases T c due to the increase in J. But 
in the overdoped region, the stripes begin to lose their 
ID character. This drives T c down since their quasi-lD 
character was the reason for the high pairing scale. 

To conclude, in the interacting 1DEG, the electron- 
phonon interaction can cause a strongly divergent super- 
conducting susceptibility with properties that are dra- 
matically different from a Fermi liquid superconduc- 
tor. Using accurate analytic techniques, we have stud- 
ied microscopic models of quasi-lD electrons coupled to 
phonons, and pointed out qualitative similarities to the 
high-temperature superconductors. These similarities in- 
clude the doping dependence of T c , the doping depen- 
dence of the superconducting pairing energy, the doping 
dependence of the isotope effect on T c , and the sign of 
the isotope effect on the spin gap. 
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